arXiv: 1509.06248v2 [hep-ph] 11 Jan 2016 


Leptoquark effects on b —> suiy and B ^ Kl^l decay processes 

Suchismita Sahoo and Rukmani Mohanta 
School of physics, University of Hyderabad, Hyderabad - 5000f6, India 

Abstract 

We study the rare semileptonic decays of B mesons induced by 6 —)• sun as well as 6 —)• sHl~ tran¬ 
sitions in the scalar leptoquark model where the leptoquarks transform as (3, 2, 7/6) and (3, 2,1/6) 
under the standard model gauge group. The leptoquark parameter space is constrained using the 
most recent experimental results on Br(i?s —)• ii~) and Br(i?rf —>■ y~) processes. Consider¬ 

ing only the baryon number conserving leptoquark interactions, we estimate the branching ratios 
for the exclusive B —)■ K^Ui/p and inclusive B —>■ Xguu decay processes by using the constraint 
parameters. We also obtain the low recoil (large lepton invariant mass, i.e., ~ mf) predictions 

for the angular distribution of i? —)• KBI~ process and several other observables including the flat 
term and lepton flavour non-universality factor in this model. 

PACS numbers: 13.20.He, 12.60.-i, 14.80.Sv 
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I. INTRODUCTION 


It is well-known that the study of B physics plays an important role to critically test the 
standard model (SM) predictions and to look for possible signature of new physics beyond 
it. In particular, the rare decays of B mesons which are mediated by flavour changing 
neutral current (FCNC) transitions are well-suited for searching the effects of possible new 
interactions beyond the SM. This is due to the fact that the FCNC transitions 6 —)■ s, d are 
highly suppressed in the SM as they occur only at one-loop level and hence, they are very 
sensitive to new physics. Recently the decay modes B —)■ which are mediated by 

the quark level transition b —)■ have attracted a lot of attention, as several anomalies 

flfl 

at the level of few sigma are observed in the LHCb experiment [iHSj. Furthermore, the 
deviation in the ratio of rates of R —)■ Kfifi over B —)■ Kee {Rk) is a hint of violation of 
lepton universality jd]. This in turn requires the careful analyses of the angular observables 
for these processes both in the low and high q‘^ regime. 

Recently various B physics experiments such as BaBar, Belle, CDF and LHCb have 
provided data on the angular distributions of R —)■ K*Rl~ and R —)■ KRl~ decay processes 
both in the low and the large recoil region except the intermediate region around ~ 
and The intermediate region is dominated by the pronounced charmonium resonance 
background induced by the decays R —)■ K{cc) —)■ KRl~, where cc = Using QCD 

factorization method the physical observables in the high recoil region can be calculated and 
the angular distribution of R —)■ Kl^l~ at low recoil can be computed using simultaneous 
heavy quark effective theory and operator product expansions in 1/Q, with Q = {mi,, \/^) 
i.e. is of the order of the 6-quark mass js, 6]. In this work, we are interested to 

study the decay process R —)■ Kll in the region of low hadronic recoil i.e. above the xjj' 
peak in the scalar leptoquark (LQ) model. We have studied the R —)■ in the large 

recoil limit in Ref. and found that the various anomalies associated with the isospin 
asymmetry parameter and the lepton flavour non-universality factor {Rk) for this process 
can be explained in this model. 

Similarly the rare semileptonic decays of R mesons with z/P pair in the hnal state, i.e., 
R ^ uu are also signihcantly suppressed in the SM and their long distance contribu¬ 
tions are generally subleading. These decays are theoretically very clean due to the absence 
of photonic penguin contributions and strong suppression of light quarks. The experimental 
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measurement of the inclusive decay rate probably be un-achievable due to the missing neu¬ 
trinos, however, the exclusive channels like B —)■ K*vv and B ^ Kvv are more promising 
as far as the measurement of branching ratios and other related observables are concerned. 
Theoretically, study of these decays requires calculation of relevant form factors by non- 
perturbative methods. 

In recent times, there are many interesting papers which are contemplated to e^lain the 
anomalies associated with the b —)■ processes, observed at LHCb experiment [l|-y , both 

in the context of various new physics models as well as in model independent ways jsl-ll 1| . 
In this paper, we intend to study the effect of scalar leptoquarks transform as (3, 2, 7/6) and 
(3, 2,1/6) under the standard model gauge group, on the branching ratio as well as on other 
asymmetry parameters in the low-recoil region of i? —)■ Kl^l~ process. We also consider 
the processes B —)■ and B —)■ Xsuu involving the quark level transitions b —)■ sui' in 

the full physical regime. It is well-known that leptoquarks are scalar or vector color triplet 
bosonic particles which make leptons couple directly to quarks and vice versa and carry both 
lepton as well as baryon quantum numbers and fractional electric charge. Leptoquarks can 
be included in the low energy theory as a relic of a more fundamental theory at some high 


energy scale in the extended SM 


12| . such as grand unihed theories 
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131, Pati-Salam 


models, models of extended technicolor M and composite models jl^. Leptoquarks are 
classihed by their fermion number [F = 3B + L), spin and charge. Usually they have a mass 
near the unihcation scale to avoid rapid proton decay, even so leptoquarks may exist at a 
mass accessible to present collider, if baryon and lepton numbers would conserve separately. 
The leptoquark properties and the additional new physics contribution to the SM have been 
very well studied in the literature QHfl- 


The plan of the paper is follows. In section II we present the effective Hamiltonian 
responsible for b —)■ processes. We also discuss the new physics contributions due to 

the exchange of scalar leptoquarks. In section III we discuss the constraints on leptoquark 
parameter space by using the recently measured branching ratios of the rare decay modes 
Bs —)■ ii~ and Bd —>■ Xs/i+p”. The branching ratio, the flat term and the lepton non¬ 
universality factor {Rk) for the decay mode B where Z = e,/i,r at low recoil 

limit are computed in section IV. In section V we work out the branching ratio oi B KvV 
process in the full kinematically accessible physical region. The branching ratio, polarization 
and other asymmetries in i? —)• K*vv process have been computed in section VI. The 
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inclusive decay process B —)■ XsW is discussed in section VII and section VIII contains the 
summary and conclusion. 


II. THE EFFECTIVE HAMILTONIAN FOR b sBr PROCESS 


The effective Hamiltonian describing the processes induced by the FCNC b 


transitions is given by 


2l| 


^e// 


4Gf 

'Tf 


■VttV, 


ts 


2 = 7 , 9,10 


2=1 


( 1 ) 


which consists of the tree level current-current operators (Oig), QCD penguin operators 
(Os-e) alongwith the magnetic and semileptonic electroweak penguin operators Og lo- 
The magnetic and electroweak penguin operators can be expressed as 


= j^(^saf,„{msPL{R) + mbPR(^L))b)F'^'" 

h') 


= ^{srPLiR)b)M , = ^{srPL(RMH7,l) . 


( 2 ) 


It should be noted that the primed operators are absent in the SM. The values of Wilson 
coefficients which are evaluated in the next-to-next leading order at the renor¬ 

malization scale /i = rrih are taken from 22|]. Here Vgqi denotes the CKM matrix element, 
Gp is the Fermi constant, a is the hne-structure constant and Pl,r = (1 T 75)/2 are the 
chiral projectors. Due to the negligible contribution of the CKM-suppressed factor VubV*^, 
there is no CP violation in the decay amplitude in the SM. These processes will receive 
additional contributions due to the exchange of scalar leptoquarks. In particular there will 
be new contributions to the electroweak penguin operators Og and Oio as well their right- 
handed counterparts Og and 0[q. In the following subsection we will present these additional 
contributions to the SM effective Hamiltonian due to the exchange of such leptoquarks. 


A. Scalar LQ Contributions to 6 —)• sl~^l effective Hamiltonian 


There are ten different types of leptoquarks under the SU{3) x SU{2) x f/(l) gauge group 


23l |. half of them have scalar nature and other halves have vector nature under the Lorentz 


transformation. The scalar leptoquarks have spin zero and could potentially contribute to 
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the quark level transition b sl^l . Here we would like to consider the minimal renor- 


malizable scalar leptoquark model 


ITj . containing one single additional representation of 


SU{3) X SU{2) X U{1), which does not allow proton decay. There are only two such models 
with representations under the SM gauge group as = (3, 2, 7/6) and Ad/^) = (^ 3 ^ 2,1/6) 


17|, which have sizeable Yukawa couplings to matter helds. These scalar leptoquarks do 


not have baryon number violation in the perturbation theory and could be light enough to 
be accessible in accelerator searches. The interaction Lagrangian of the scalar leptoquark 
AT/ 6 ) -vv-ith the fermion bilinear is given as jl^ 


=gRQLA^^^^^lR + h.c., 


( 3 ) 


where Ql is the left handed quark doublet and Ir is the right-handed charged lepton singlet. 
After performing the Fierz transformation and comparing with the SM effective Hamiltonian 
([I]), one can obtain the new Wilson coefficients as discussed in Ref. [l^ 

{9R)sl{gR)li 


^NP _ f _ 

'-'9 ~ '-'10 ~ 


■yAfP 


vr 


( 4 ) 


2V2GfaVtbV* 

Similarly, the Lagrangian for the coupling of scalar leptoquark Ad/^) to the SM fermions is 
given by 

£(1/6) = ^^4Ad/6)t£ + h.c., with A = (5) 

where T 2 is the Pauli matrix and consists of operators with right-handed quark currents. 
Proceeding like the previous case one can obtain the new Wilson coefficients as 

TT {gL)sl{gL)ti 


/~iiNP _ _ 

'-^9 ~ '-'10 ~ 


'i/NP 


( 6 ) 


2^G,aVaVi ’ 

which are associated with the right-handed semileptonic electroweak penguin operators Og 


and O', 


10 - 


III. CONSTRAINT ON THE LQ PARAMETERS 


After having the idea of possible scalar leptoquark contributions to the b —)■ sll processes 
we now proceed to constraint the LQ couplings using the theoretical 2^ and experimental 


branching ratio 


25 


27| of t process. This process is mediated by 6 —)■ s/i/i 


transition and hence well-suited for constraining the LQ parameter space. In the SM the 
branching ratio for this process depends only on the Wilson coefficient Cig. However, in 
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the scalar LQ model there will be additional contributions due to the leptoquark exchange 
which are characterized by the new Wilson coefficients and C'iq^ de pen ding on the 
nature of the LQs. Thus, in this model the branching ratio has the form 71. Il9| 


Bt{Bs 


* 12 

^l\ ’ m’ ts\ 


/~iSM 1 /~i 

Oin O 


NP 


c 


NP 


4m?, 

_P 

Ml. 


m 


which can be expressed as 

Br(5, ^ /r+/i-) = Br^^ 


1 + 


r<NP r"NP 
'^10 *-^10 


CSM 

*-^10 




1 + 


where Br"^^ is the SM branching ratio and we dehne the parameters r and (f)^^ as 


re 


r^NP r^'NP 
_ *^10 *-"10 
“ rSM 
*-"10 


Now comparing the SM theoretical prediction of Br(i?s —)■ /i/i) 24 1 

Br(5, ^ /r+/i“)lsM = (3.65 ± 0.23) X 10“^ 
with the corresponding experimental value 

Br(S, ^ /i+/i”) = (2.9 ± 0.7) x 10-^ 


( 8 ) 


( 9 ) 


( 10 ) 


( 11 ) 


one can obtain the constraint on the new physics parameters r and . The constraint on 


the leptoquark parameter space has been extracted in jj, ll9|] from this process, therefore, 
here we will simply quote the results. The allowed parameter space in r — plane which 
is compatible with the la range of the experimental data is 0 < r < 0.1 for the entire range 
, i.e., 

0 < r < 0.1 , for 0 < < 27i . (12) 


of p^^ 


However, in this analysis we will use relatively mild constraint, consistent with both mea¬ 
surement of Br(i?s as 

0 < r < 0.35 , with n 12 < p^^ < 37r/2 . (13) 

It should be noted that the use of this limited range of CP phase, i.e., (7r/2 < p^^ < 37r/2) 
is an assumption to have a relatively larger value of r. These bounds can be translated to 
obtain the bounds for the leptoquark couplings as 

0 < < 5 X 10-9 GeV-2 for 7r/2 < p^^ < 37r/2 . (14) 
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After obtaining the bounds on leptoquark couplings, we now proceed to study the decay 
processes B —)■ Kll and B —)■ K^*\Xs)uP and the associated observables in the following 
sections. 


IV. B Kl+l- PROCESS IN THE LOW-RECOIL LIMIT 


The transition amplitude for the B —)■ Kl^l~ decay process can be obtained using the 
effective Hamiltonian presented in Eq. ([T]). The matrix elements of the various hadronic 
currents between the initial B meson and the final K meson can be parameterized in terms 


of the form factors /o, fr and /+ as 


28| 


{K (k) \s-f>^b\B (p)) = /+ (g2) (p + kf + [fo (q^) - /+ (g^)] 
{K (k) \sa^^b\B (p)) = i ^ [(p -)_ q'^ _ ('p _j_ ^ 


m. 


m 


I£ tj. 


q" 




(15) 


(16) 


rriB + rriK 

where p, k are the four-momentum of the H-meson and Kaon respectively and q = p — k 
is the four-momentum transferred to the dilepton system. Furthermore, using the QCD 
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operator identity 


an improved Isgur-Wise relation between and /+ can be obtained as 


id"" {siaf,„b) = -rrib (s'y^.b) + id^ (sb) - 2 (si^^b^ , 


/t (g^/i) = 


ms {mB + t hk) 

q2 


-K 


(a)/+(9")+o(AT 

\mbj 


(17) 


(18) 


where strange quark mass has been ne glec tec 
B —>■ Kl^l~ process in low recoil limit 28 


Thus, one can obtain the amplitude for the 




TT 


3lj, after applying form factor relation fflSj) as 

(19) 


A(B^ Kl*r) = i^V,tV:jAf)[Fv1^{h^l) + Fap“(IwI) + Fp(h!,l)], 


where 


Fa=C- 


tot 
10 5 


Fv = cr + K 


2mbmB 

q2 




Fp = -mi 


1 + 


m 


B 


m 


K 


fo 

u 


( 20 ) 


y-itot 
^10 • 


In Eqn. ([20]), and , where and 

are the new contributions to the Wilson coefficients arising due to the exchange of lep- 
toquarks and the effective Wilson coefficients are given in Ref. 321. The corresponding 
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differential decay distributions is given by 


d^Ti [B KI+1-] 


= ai [q^) + Cl (g^) cos^ Oi , 


( 21 ) 


dq^d cos 9i 

where 9i is the angle between the directions of B meson and the in the dilepton rest 
frame. The expressions for the q^ dependent parameters a^, q are presented in Appendix A. 
Thus, the decay rate for the process B —)■ Kl^l~ can be written as 


f*q^ 

^max 


T, = 2 


■ 

^rmn 


dq^ 


0^1 + -Cl 


( 22 ) 


Another useful observable known as the flat term is defined as 

2 


' ^max 


F^h = 


^min 


dq‘ (a, + ci), 


(23) 


where the hadronic uncertainties are reduced due to cancellation between the numerator 
and denominator. It should be noted that the lepton mass suppression of (a/ + q) follows 
as oc mf, hence, it vanishes in the limit mi —)■ 0. 

After obtaining the expressions for branching ratio and the observable Fi, , we now proceed 
for numerical estimation for B —)■ Kl^l~ process in the low recoil region. In our analysis we 
use the following parametrization for the dependence of form factors /j (z = +,T, 0) as 


28|,l33| 


fi (s) = 


/.(O) 


res,I 


1 + b\ (^z (s) -z{0) + ^{z {sf - z (0)^)^ 


1 — s/m\ 

where we have used the notation = s. The z(s) functions are given as 
y/r+ - s - ^T+ - To 


(24) 


z(s) = 


\/t-+ - s + \/^+ - To ’ 


To = , 


r± = (ms ± mxY 


The values of fi{0) and b\ are taken from 


For numerical evaluation, we have used the particle masses and the lifetimes of B meson 
from 3^. For the CKM matrix elements, we have used the Wolfenstein parametrization 


with values A = 0.814+°;°23, A = 0.22537 ±0.00061, p = 0.117 ±0.021 and fj = 0.353 ±0.013 
and the hne structure coupling constant a = 1/137. With these input parameters, the 
differential branching ratios for B^ —)• K^e^e~ (left panel), B'^ —)■ (right panel) 

and —)■ (lower panel) processes with respect to high g^, both in the SM and in 

the leptoquark model are shown in Fig. 1 for leptoquark and in Fig. 2 for 

The grey bands in these plots correspond to the uncertainties arising in the SM due to the 

















uncertainties associated with the CKM matrix elements and the hadronic form factors. The 
green bands correspond to the LQ contributions. For 5 —)■ process, we vary the values 
of the leptoquark couplings as given in Eq. (14) and for B —)■ Kee and B Ktt processes 
we use the limits on the LQ couplings extracted from B^ —t Xae'^e~ and Bg —)■ r+r“ 


0 


processes |Y1J as 


and 


0 < I fee I < ^ ^ ^Q- 

- Ml - 


GeV 


-2 


g ^ \{SR).A!}R)kr\ < J 2 ^ jg-S 

Ml 


(25) 


(26) 


Since the leptoquark couplings are more tightly constrained in 6 ^ s/i/i transitions, the 
deviations of the branching ratios in the LQ model from the corresponding SM values are 
found to be small. For B —)■ Kee and B —)■ Ktt these deviations are found to be signihcantly 
large. The bin-wise experimental values are shown in black in 5 —)■ Kfifi process. From 
these hgures it can be seen that the observed experimental data can be explained in the 
scalar LQ model but the deviation from the SM branching ratios are more in the A^/®) 
model. For the other observables in 5 —)■ Kll processes we will show the results only for 
leptoquark model. In Fig. 3, we have shown the lepton non-universality factors 
(left panel) (he. the ratio of branching ratios of 5 —)■ KfK and B —)■ .^e’''e“), (right 
panel) and (lower panel) variation with high From the hgure one can see that there is 
signihcant deviations in the lepton-flavour non universality factor from their corresponding 
SM values in all the above three cases. The flat term for the B^ —)■ K^R ^(left panel) 
and B^ —)■ K^t^t~ (right panel) decay processes in the low recoil region are presented in 
Fig. 4 for In this case there is practically no deviation in B ^ Kfifi whereas there 

is signihcant deviation in i? Ktt process. The integrated branching ratios, hat terms 
and the lepton havour non-universality factors for the B —)■ Kll processes over the range 
G [14.18,22.84] are given in Table 1. The hat term for B —)■ Ke'^e~ process has been 
found to be negligibly small {Ffj ~ 0(10“^)) due to tiny electron mass. In the low recoil 
region, the process having tan lepton in the hnal state has signihcant deviation from the 
SM. 

The integrated branching ratio for B^ —)■ Kfui process in the range G [15,22] GeV^ 
has been measured by the LHGb Gollaboration Ij and is given as 

Br(5° ^ K^fifi) = (6.7 ± 1.1 ± 0.4) x 10"® . (27) 
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Our predicted value in this range of is found to be 


(SM) 

LQ model) 
LQ model) 


Br(5° ^ iLVAi) = (8-35 ± 0.5) x 10“®, 

= (8.34- 9.26) X 10"® . 
= (8.34-15.6) X 10"® . 


(28) 


The predicted values of the branching ratios are slightly higher than the central measured 
value but consistent with its l-cx range. 
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FIG. 1: The variation of branching ratio for —?■ K^e^e~ (left panel), B —>• (right panel) 

and B —?■ (bottom panel) with respect to high for AT/6 ) lQ. The grey bands correspond 

to the uncertainties arising in the SM. The q^-averaged (bin-wise) 1 — a experimental results for 
B —)• K/u/j, process are shown by black plots, where horizontal (vertical) line denotes the bin width 
(1 — cr error). 
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FIG. 2: Same as Fig-1 for LQ exhange. 

V. B ^ Kvu PROCESS 


The B ^ Kvv process is mediated by the quark level transition b —)■ suiy and the effective 
Hamiltonian describing such transition is given as 351 

«.// = + ^rOr) + h.c., (29) 

where 

2 2 

{sip^PLh) (1 - 75) J^), {sl^,PRb) (h7^ (1 - 75) v) , (30) 

are the dimension-six operators and C'[ ^ are their corresponding Wilson coefficients. The 
coefficient C'^ has negligible value within the standard model while C'[ can be calculated by 
using the loop function and is given by 

Cl = -X{xt)/sm^e^. (31) 
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FIG. 3: The variation of lepton non-universality (left panel), (right panel) and (bottom 

panel) in low recoil region due to LQ exchange. 


The necessary loop functions are presented in Appendix B. Th e decay distribution with 
respect to the di-neutrino invariant mass can be expressed as 


dr 


ds 


B 


Gy 

2567r5 




V |2 


(32) 


where rfij = mi/rriB and sb = s/m?B. The decay rate has been multiplied with an extra 
factor 3 due to the sum over all neutrino flavours. It should be noted that in Eq. fl3^ 
is the new Wilson coefficient arises due to the exchange of the leptoquark Ad/®). In order 
to hnd out its value, we consider the new contribution to the effective Hamiltonian due to 
the exchange of such leptoquark which is given as 


A /fi'l 


(33) 
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FIG. 4: The variation of flat term for —)• (left panel) and —)• (right panel) 

with high for LQ. 


Oservables SM predictions Values in LQ model Values in LQ model 


Br(B0 ^ K^e+e-) 

(1.005 ±0.06) 

(1.004 - 1.5) 

(1.005 - 1.88) 

Br(B° ^ ATV+d") 

(1.01 ±0.06) 

(1.01 - 1.12) 

(1.008 - 1.89) 

Bt{B^ K^t+t-) 

(1.21 ±0.73) 

(0.99 - 2.07) 

(1.2 -4.2) 

{K) 

1.0035 

0.75 - 1.00 

1.0035 

{R^k^) 

1.21 

0.98-1.85 

(1.2 - 2.3) 

(RW) 

1.198 

0.98-1.85 

(1.2 - 2.2) 

{Fh) 

1.75 X 10"^ 

(1.74 - 1.75) X 10-'^ 

(1.73 - 1.75) X 10' 

{F^h) 

7.5 X 10“^ 

(7.4 - 7.55) X 10-3 

(7.4 - 7.5) X 10-' 

(Fh) 

0.89 

0.8 - 1.38 

0.88-0.89 


TABLE 1: The predicted values for the integrated branching ratios (in units of 10“'^), flat terms and 
lepton non-universality factors in the range G [14.18,22.84] GeV^ for the B —?■ Kl^l~ process, 
I = e, fj,, T. 


Comparing Eqs. fl29|) and fl33|) . one can obtain the new Wilson coefficient as 

QU \ _ __ {9 l) su{9 l) bu 

“ 2V2GFaVtbV,l Ml,,,,, 


(34) 
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For numerical estimation, we use the B ^ K form factor evaluated in the light cone 

n 

sum rule (LCSR) approach [371] as 


fn<in = 


n 


+ 


r2 


(35) 


1 — q‘^/m‘1 (1 — /m\Y ’ 

which is valid in the full physical region. Furthermore, in contrast to i? ^ Kl^l~ process, 
which has dominant charmonium resonance background from B —)■ K^J/tjj) —)■ Kl^l~, 
there are no such analogous long-distance QCD contributions in this case as there are no 
intermediate states which can decay into two neutrinos. For the b —)■ suu LQ couplings 
we use the values as we used for b —)■ s/i/i as these two processes are related by SU{2)l 
symmetry. The variation of branching ratio with respect to ss in the full physical regime 
0 ^ Ss ^ (1 — friKY is shown in Fig. 5 and the predicted branching ratio is given in Table 
II, which is well below the present upper limit Br(i?2 —)• Kvv) < 4.9 x 10“® 34|. 



FIG. 5: The variation of branching ratio ol B ^ Kuv with respect to the normalized invariant 
masses squared sb in the SM and Ad/6) j^odel. The grey band corresponds to the uncertainties 
arising in the SM. 


VI. B K*Tyv PROCESS 

The study of i? —?■ K*u9 is also quite important as this process is related to B ^ K*^fi 
process by SU{2)l and therefore, the recent LHCb anomalies in i? —)■ K*ij,ij, would in 
principle also show up in i? —)■ K*i>v. The experimental information about this exclusive 
decay process can be described by the double differential decay distribution. In order to 
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compute the decay rate, we must have the idea about the matrix element of the effective 
Hamiltonian fl2^ between the initial B meson and the hnal particles. Due to the non¬ 
detection of the two neutrinos, experimentally we can’t distinguish between the transverse 
polarization, so the decay rate will be the addition of both longitudinal and transverse 
po 


35 


arizations. The double differential decay rate with respect to and cos^ is given by 

(PV 


3 dV'J' .Or, 3 dr L 2 

sm^ 9+ 


dssd cos 9 AdsB 2dsB 

where the longitudinal and transverse decay rate are 


cos 9, 


dT, 


= 3m||HoP, 


dTn 


= 3m^^(|HxP+|H||p). 


(36) 


(37) 


dsB dsB 

The transversality amplitudes ^±,||,o in terms of the form factors and Wilson coefficients are 
listed in Appendix C. 

The fractions of K* longitudinal and transverse polarizations are given as 

d^B^r/ dsB 


and the K* polarization factor is 


Fl,t — 


OiK* — 


dV/dsB ’ 


2Fl 

Ft 


1 . 


The transverse asymmetry parameters are given as 39|, | 

„ _ -2Re{AyAl) ^ 

^ l^iP+IAlP’ l^iP+IAlP' 


( 38 ) 


( 39 ) 


( 40 ) 


However, one can’t extract from the full angular distribution of i? —)■ as it is not 

invariant under the symmetry of the distribution function and requires measurement of the 
neutrino polarization. So it can’t be measured experimentally at B factories or in LHCb. 
The transverse asymmetry Arp' is theoretically clean and could be measurable in Belle-H. 
For numerical evaluation we use the dependence of the B —)• K* form factors 


l/(g^), Ai{q 2 ), A 2 {q‘^) from |^. The variation of the branching ratio of 5 —)■ K*vv 


with respect to the neutrino invariant mass, ss is shown in Fig. 6. Fig. 7 contains the 
longitudinal and transverse polarizations of K* verses sb- The polarization factor and the 
transverse asymmetry variation with respect to ss in the full region are shown in Fig. 8. 
Although there is certain deviation found between the SM and LQ model predictions for 
the branching fraction, but no such noticeable deviations found between the SM and LQ 
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predictions for the longitudinal/transverse polarizations, transverse asymmetry parameters 


( 0 \ 

A>p . The integrated valnes of branching ratio over the range sb ^ [0, 0.69] are presented 
Table II, which are well below the the present npper limit < 5.5 x 10 

34|. 


m 


-5 





FIG. 6; The variation of branching ratio of i? —)• K*vv with respect to the ss in the SM and A^/®) 
LQ model. The grey band corresponds to the uncertainties arising in the SM. 




Sb Sb 


FIG. 7: The variation of longitudinal (left panel) and transverse (right panel) polarization of K* 


with SB- 
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Sb 


Sb 


FIG. 8: The variation of K* polarization factor (left panel) and the transverse asymmetry (right 
panel) with respect to ss. 


VII. B Xsvv 


The inclusive decay B —>■ Xgvv is dominated by the Z-exchange and can be searched 
through the large missing energy associated with the two neutrinos. This decay mode is 
theoretically very clean, since both the perturbative and the non-perturbative corrections 
are small. So these decays do not suffer from the form factor uncertainties and thus, are 
very sensitive to the search for new physics beyond the SM. The decay distribution with 
respect to Sf, = s/ml can be written as 


dr 

dsh 


a^G 




X 


1287r5 
3sb(l + ml 


Sb - ^ms 


+ A(l,m,,Sb) 


(41) 


'\Cl? + \C^R?' 

where rhj = m^/mb and fi:(0) = 0.83 is the QCD correction to the b —)■ suu matrix element 
431]. The full kinematically accessible physical region is 0 < s;, < (1 — rfis)'^- In Fig. 9, 
we have shown the variation of the branching ratio with respect to Sb and the integrated 
branching ratio values over the range Sb G [0, 0.96] both in the SM and in the leptoquark 
model are presented in Table II. 


We dehne the ratio of branching ratios as 


36|, 


—)■ Ki'h') 
Bt{B Xsuu) ’ 


(42) 
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and 


Rk* = 


Bt{B —)■ K*uv) 


(43) 


Br(i? —)■ Xgiyu) 

and the variation of Rk and Rk* with respect to sb in the fnll kinematically allowed region 


is shown in Fig. 10. In this case also no deviation fonnd between the SM and leptoqnark 
predictions. 



FIG. 9: The variation of branching ratio of B ^ XsW with respect to the Sf,. 




Sb Sg 


FIG. 10: The variation of Rk (left panel) and Rk* (right panel) with respect to sb- 


VIII. CONCLUSION 

In this paper we have stndied the effect of scalar leptoqnarks on the rare semileptonic 
decays of B meson. In particniar, we focns on the decay processes B —)■ Kl^l~ in low recoil 
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Observables 

SM prediction 

Values in AB/®) LQ model 

Br(R0 ^ 

(4.9 ±0.29) X 10"® 

(3.6 - 5.2) X 10-6 

Br(R0 ^ K*vu) 

(9.54 ±0.57) X 10-® 

(7.02 - 10.13) X 10-6 

Br(i? —)• Xgvy) 

(2.98 ±0.18) X 10-5 

(2.2 - 3.17) X 10-6 

Rk 

0.164 

(0.163 - 0.164) 

Rk* 

0.32 

0.32 


TABLE II: The predicted branching ratios for B —?• {K, K*, Xs)i'i' processes and Rk,k* for B —>• 
Xgvv in their respective full physical ranges. 

limit and the di-neutrino decay channels B —)■ K^*\Xs)i'v. The leptoquark parameter space 
is constrained by considering the recently measured branching ratios of —)■ and 

Bd —)■ Xs/i+/r“ processes. Using the allowed parameter space we predicted the branching 
ratio, lepton non-universality factors and the flat terms for the B —)■ Kl^l~ process in the 
low recoil region. We found that the measured branching ratio can be accommodated in 
the scalar leptoquark model. We have also calculated the branching ratios of 5 —)■ 
and B —)• Wz/P processes. The predicted branching ratios for B —)■ processes are 

well below the present upper limits. The polarization of K* and transverse asymmetry for 
B —)■ K*i>v are also computed using the constraint leptoquark parameters. However, we 
found no deviation between the SM prediction and the LQ results for different polarization 
variables and the transverse asymmetry parameter. 
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Appendix A: a; and q functions in B —> Kll process 

The ai and q parameters in the decay distribution of the B —)■ Kl^l~ processes can 
be expressed as 

^ ^2 = I -mK + (f) Re {FpFD 

+ Aw?imWFA\^ + q^\Fpf‘, (Al) 
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with 


and 


Q 


^o^/AA/^ 


=-^fUlFAl^ + lFyn , 





\ 4, 4|4ot22| 22, 22\ 

X = m^ + mj^ + q - 2 [m^rrij^ + m^q +m^q ). 


(A 2 ) 


(A3) 


Appendix B: Loop functions 


The loop function X(xf) in Eq. (ET]), including correction 0{as) at the next-to-leading 
order in QCD, is given by j44J . l45 | 


a.. 


X{xt) = Xo(xt) + —Xi(xt), 

47r 


where 


Xo(xt) = 


Xt 


2 + Xt 3xt-Q 

+ 7^-77 In Xt 


1-xt (1 - xty 


(Bl) 


(B2) 


and 


Xi{xt) = 

+ 


29xt — Xt — 4xt Xt + 9xt — xf — xf 


3(1(l-a;t)3 

8xt + 4:Xt + Xt — Xt , 2 


f '^f '^f 1 

\xixt 


In Xt — 


Axt — Xt 


L2{1 - Xt)+ 8xt^^^^^^\nxt,. (B3) 


2(1 -xt)^ (l-Xi)^""'* dxt 

In Eqs. flBlllB3p . the parameters used are defined as Xt = m^/rri^, = ^‘^/rri^ with 
fi = 0{mt) and L 2 (l - Xt) = 


Appendix C: Transversity amplitudes for B —)• K*vv process 


The transversality amplitudes A_lj|^o for B K*vv process are given as 

A±(ss) = 2N \/2A^^^(l, sb){C'1 + n ~ ^ i ’ 

(1 +mK*) 

A||(sb) = —2NV2{1 + mK*){C'[ — C'y)Ai{sB), 

N{Cl - Cy) 


^o(ss) — 


mK*\fSB 


(1 - rriK* - Si?)(l + mK*)Ai{sB) - X{l,mj^*,SB) 


A2{sb) 
1 + friK* 


(Cl) 

(C2) 

(C3) 
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with 


N = VtbV, 


ts 


3 • 2i07r5 


Sb 




1/2 


(C4) 


The various form factors V{sb), Ai(sb), ^ 2 (^ 5 ) associated with B —)■ K* transition in Eqs. 
(C1-C3) are dehned as 

{K* (pk*) \s'y^PL,Rb\B (p)) = ie^^ai3e''*p°q^ _ — T ^ | {ruB + mK*)elAi{sB) 


rriB + rriK* 


-{e*-q){2p- q)^ 


^2(55) 


2mK* 


mB + rriK* s 

where q = pi+ + pi- and is the polarization vector of K*. 


(e*-?) [ 413 ( 55 )-Alo(sB)]g/. , (C5) 
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